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I. First overtone bands from a disk or from a wind? 
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Abstract. We observed the near infrared emission in the 
wavelength range 2.28-2.5 /im from the peculiar B[e]-star 
MWC 349. The spectra contain besides the strong IR contin- 
uum the first overtone CO bands and most of the hydrogen re- 
combination lines of the Pfund series, both in emission. We 
also modeled the spectra. The Pfund Unes have a gaussian pro- 
file with a FWHM of ^ lOOkm/s, and it turned out that their 
emission is in LTE and optically thin. To explain the CO bands, 
several scenarios were investigated. We found that the CO band 
heads are formed under LTE and that the gas must have a tem- 
perature of 3500 to 4000 K. The width of the 2 band head 
indicates kinematical broadening of 50 to 60km/s. We can ob- 
tain fits to the measured spectra assuming that the CO gas has a 
column density of 5 • 10^'' cm^^ and is located either at the in- 
ner edge of the rotating circumstellar disk. In this case, the disk 
must have a bulge which partly blocks the radiation so that the 
observer sees only a sector on the far side where the radial ve- 
locities are small. Or the CO emission originates in a wind with 
gaussian line profiles. Both fits are of equal quality and satis- 
factory. In a third alternative where the fit is less convincing, 
the CO emission is optically thin and comes from an extended 
Keplerian disk. 

Key words: Line: formation - circumstellar matter - Stars: 
emission-line. Be - Infrared: stars - Stars: individual: 
MWC 349 



1. Introduction 

After the discovery of CO band head emission in the BN object 



by Scoville et al. (1979) a large number of further detections in 



other y oung st ellar object s (YSO) follow ed (se e e.g. Geballe & 



1996;Najitaetal. 1996) 



Persson 1987 Carr 1989 Chandler et al. 1993 ; Greene & Lada 



Several scenarios were discussed to explain the origin of 



the hot (2500-5000 K) and dense (n > 10^^ cm 



nent (for some examples see Calvet et al. 1991 



^) CO compo 
Martin 



1997), 
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but the most likely location is a neutral disk or wind (Carr 1989 



Chandler, Carlstrom, & Scoville 1995). The disk model is espe 



cially supported by high resolution spectroscopic observations 
of the CO 2^0 band head. The shape of this band head shows 
for several YSO the kinematic signature of Ke plerian rotation 
(Carr et al. |1993| ; Carr |l995| ; Najita et al. |1996| ) and is a power- 
ful tracer for the existence of a circumstellar disk around young 
stellar objects. 

Also the peculiar B[e]-star MWC 349 shows the first over- 
tone CO bands in emission. These bands were observed first by 
Geballe & Persson (1987) with a velocity resolution of about 
460km/s. 

MWC 349 is a binary system consisting of the main com- 
ponent MWC 349 A, classified as a B[e]-star, and the BO III star 
companion, MWC 349B, localized 2" west of MWC 349A. We 
are only interested in the main component, in the following ref- 
ered to as MWC 349. Its evolutionary state is still unclear. It 
shows some characteristics of a pre-main sequence B[e]-type 
star as well as characteristics of a B[e] supergiant (e.g. Lamers 
et al. 



1998) 



Cohen et al. ( |1985[ ) determined its distance (1.2 kpc), bolo- 
metric luminosity (^ 3 • IC* Lq) and the visual extinction to- 
wards it (Ay ~ 10 mag of which 2 mag might be circumstel- 
lar). 

The existence of a bulge of circumstellar dust around 



MWC 349 has been known since long ago (e.g. Geisel 1970) 



The proposition that MWC 349 also has a disk was supported 
by obs ervations of double-peaked emission lines (Hamann & 
Simon 1986 , 1988 1 and b y IR speckle interferometry (Leinert 
Mario tti et al. 



1986 



1983) which revealed a disk-like structure 



of the dust emission in the east-west direction, seen nearly edge 
on. An accumulation of neutral gas and dust in the equatorial 
plane of the star is also expected to be responsible for the bipo- 
lar structure of the optical ly thic k wind zone seen in the VLA- 
map of White & Becker (1985). The mass loss rate found for 
a 50 km/s wind velocity is ^ 1.2 • 10^'^ MQyr^^ (Cohen et al. 



1985). 



Another indicator for the circumstellar disk are the strong 
hydrogen recombination maser lines in the mm and submm 
range which also show the characteristic double-peaked pro- 
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files (e.g. Martm-Pintado et al. 1989). Recombination lines at 
different wavelengths sample different regions. With decreas- 
ing quantum number n, i.e. increasing frequency, one sees ion- 
ized gas closer to the star The fact that the rotational velocity 
of recombination lines with decreasing n displays a systemati- 
cal increase led to the assumption that the double-peaked maser 
emission comes from the ionized atmosphere of a Keplerian ro- 
tating disk around a 25-30 M q sta r (Thum, Martm-Pintado & 
Bachiller 1992; Thum et al. 1994) whose free-free and free- 



Pa„ 



bound emission was modelle d by Kraus et al. ( [2000| ). In ad- 
dition, Rodriguez & Bastian ( 1994 ) determined the inclination 
angle of the disk towards the line of sight to 15° ± 5°. 

In this paper, we present new low resolution observations of 
the total first overtone CO band emission as well as high reso- 
lution observations of the 2 — > and 3^1 band heads. We 
discuss the probable location of the hot CO gas by modeling 
the emission for several scenarios. 

2. Observations 

In 1997 we observed MWC 349 with the UKIRT telescope in 
the NIR range (1.85 — 2.5 /im) with a spectral resolution of 
^ 330km/s (Fig.[I^). The spectrum contains several prominent 
hydrogen and helium recombination lines as well as the first 
overtone CO bands, all of them in emission. 

Table 1. Identified lines in the NIR spectrum of MWC 349 
(Figs. |l^ and b). 



number 


wavelength [/im] 


line identification 


1 


1.86904 


He I 


2 


1.87561 


Paa 


3 


1.89282 


Cl, Cal 


4 


1.90867 


He I 


5 


1.94508 


Br5, Cai 


6 


1.95445 


Nl 


7 


1.97470 


Sil, Cl,Nl 


8 


1.98702 


Cal, Cl 


9 


1.99318 


Cal, Sil, Cl 


10 


2.00130 


Cl 


11 


2.05917 


He I 


12 


2.08819 


(?)SnI, Fell 


13 


2.11194 


He I 


14 


2.13831 


Mg II, C I, Sir 


15 


2.14359 


Mgll, Call, [Felll] 


16 


2.16611 


Br7 


17 


2.20687 


Na I doublet, Sil, (?)OlI 


18 


2.21741 


Cl 


19 


2.24113 


(?) Fell 



In two sessions in 1998 we reobserved MWC 349 in the 
spectral range 2.285-2.341 /im with higher spectral resolution 
(10-15km/s). The observations consist of four seperate sub- 
spectra with a small overlap at their edges. They are slightly 
offset relative to each other and at their low wavelength end, 
the fluxes are systematically overestimated. We calibrated the 
subspectra by smoothing them to the resolution of the spectrum 
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Fig. 1. For better visualization, the total NIR spectrum of 
MWC 349 observed with UKIRT with a spectral resolution of 
~ 330km/s (upper panel) is split in two parts: Part 2 (bottom 
panel) contains the first overtone CO bands (indicated by the 
arrows) and most lines from the Pfund series of the hydrogen 
atom (vertical bars), the remaining spectrum (Part 1, middle 
panel) contains several emission lines listed in Table |l|. 



in Fig. |l]c, then scaled each spectrum with a factor near ^ 1 and 
combined them to the final plot of Figure ^. 
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Unfortunately, the 3 — > 1 band head falls into the over- 
lap region of two subspectra, so its strength is rather uncertain 
(marked by the dashed line within the figure). This fact turned 
out during the modeling to be a disadvantage. 
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Fig. 2. Resulting calibrated high resolution spectrum of the first 
two CO band heads from MWC 349. The dashed regions mark 
the overlap for each two subspectra where the calibration is 
rather uncertain. 



3. Results 

3.1. Hydrogen recombination lines - the Pfund series 

As can be seen in Figs. |l|c and ^ the spectrum of MWC 349 in 
the wavelength range from 2.25 to 2.5 /xm contains beside the 
first overtone CO bands also the Pfund series (transitions of the 
form n ^ n' = 5, deno ted as Pf(n)) of the hydrogen atom. 
Hamann & Simon ( 1986 ) already identified some of them in 
their velocity-resolved infrared spectroscopy of MWC 349, but 
their resolution of ^20km/s was lower than ours, so they could 
only resolve lines up to Pf(34). Their FWHM values of several 
unsmeared lines were about lOOkm/s and agree with ours. 

Many of the other emission lines they observed are double- 
peaked which is a hint for a rotating medium. At first glance 
we cannot exclude the Pfund lines to take part of the rotation 
because they are superimposed on the CO bands, but their pro- 
file looks more than a gaussian than a double-peaked one. A 
gaussian velocity component of about 50km/s, however, is not 
a surprising feature for MWC 349. Also the cm and mid IR hy- 
drogen recombination lines as well as the 'pedestal' feature of 
the double-peaked mm recombination maser lines show such a 
velocity co mpone nt (Smith et al. 1997 ; Thum, Martm-Pintado 
& Bachiller[l99^. It is often ascribed to the wind of MWC 349. 



To model the hydrogen emission, we assume the Pfund 
lines to be optically thin. Their intensity is given by 



(1) 



To evaluate it, we need the number density iV„ of the atoms in 
level n, the Einstein coefficients, An.^, and the profile function 
of the hydrogen gas, ^ii{v). 

In view of the measured line widths, we choose a gaussian 
profile function defined as 

2' 



$h(i^) = 



1 



exp 



1^0 



(2) 



V is the most probable velocity which is in the case of the ther- 
mal velocity given by w = yj2kT jm for particles having a 
Maxwellian velocity distribution. In our case, v ~ 50km/s. We 
also take into account the local standard of rest velocity of 
- 8km/s (Thum, Martm-Pintado & Bachiller |l992t Thum et 
al. 



1995) 




2.3 



2.31 



2.32 
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Fig. 3. Spectrum of optically thin Pfund lines in LTE, with an 
isothermal electron temperature of Tc — lO'^K and a gaus- 
sian velocity component of 50km/s. The spectral resolution to 
which the lines were smoothed are 330km/s (dashed line) and 
15 km/s (solid line). At short wavelengths, the lines blend into 
a 'continuum', especially for the low resolution. 

The Einstein coefficients, A„„/, can be calculated, follow- 



ing Menzel & Pekeris ( 1935 ). If we restrict ourselves to Pfund 
lines with n > 18 (see Fig.lllc), they can be approximated by 



A 



-71,5 



3.023- 10^ n 



9 „-5.005 



(3) 



The number density, Nn, of the level n follows from Saha's 
equation 



2 (27rmefcre)3/2 



(4) 
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Fig. 4. Subtraction of the synthetic continuum and Pfund line 
emission (solid line, middle panel, pure continuum is shown 
by the dashed line) from the observations (upper panel) leads to 
the almost pure CO first overtone band emission (lower panel). 



Detailed modeling of the Pfund line spectrum performed 
by Kraus (2000) showed that the emission arises within the 
innermost parts of the ionized wind where the electron den- 
sities are highest. Within these regions it is expected that < 
Electron densities as high as 10^ cm^'^ have already 



been found by Sti-elnitski et al. ( |1996| ) and Thum et al. ( |1998[ ). 
For such high electron densities the emission is clearly in LTE, 
i.e. the 6„-factors are ^ 1. Therefore, non LTE effects are neg- 
ligible. 

In our case of optically thin Pfund lines in LTE and with a 
constant electron temperature the intensity ratios of two neigh- 
bouring Pfund lines is according to Eq. (|l|) independent of the 
electron density 



'Pf(n+1) 
Ipf(n) 



5.005 



• exp 



hRy 2n + l 

kTc {v? + nY' 



(5) 



and always smaller than 1, which means that the intensity 
of the Pfund lines decreases with increasing quantum num- 
ber. This behaviour can qualitatively be seen in Fig.^ where 
we plotted the synthetic Pfund lines which we smoothed to 
the spectral resolution of 10-15km/s (solid line) and 330km/s 
(dashed line), respectively. The wavelength difference between 
the neighbouring Pfund lines decreases with increasing n lead- 
ing to a blend of the individual lines. The onset of this blend 
depends on the width of the lines and on the spectral resolution 
and leads to a hydrogen 'continuum'. 



Table 2. Parameters for hydrogen used for the calculations of 
the Pfund lines shown in Fig. Q (middle panel). 





Tc [K] 


Ugauss [km/s] 


A'^'^'im&g] \ N'^dV [cm.-'] 




< 10" 


~ 50 


10 6.35 ■ 10"" 



The interstellar visual extinction towards MWC 349 is 
about A}^^ ~ 10 mag and has to be taken into account. Fit- 
ting the Pfund lines with the parameters displayed in Table ^ 
leads to the quantity J N^dV, where V is the volume of the 
emitting region. 

The middle panel of Figure Q contains the modelled Pfund 
series added to the continuum emission of MWC 349 (dashed 



line) taken from Kraus et al. (2000). These two components 
were subtracted from the observations (upper panel) and the 
resulting spectrum (lower panel) can mainly be ascribed to the 
CO first overtone band emission. To fit this observed CO spec- 
trum is the aim of the following sections. 



where iVj and Nc are the number densities of the protons and 
electrons, To and nic are the electron temperature and electron 
mass, gn = 2n? is the statistical weight of the state n, and Ry is 
the Rydberg constant. The 6„-factors which describe the devi- 
ation from LTE and which depend on the electron temperature 
and densi ty wer e taken from the electronic Tables of Storey & 
Hummer ( 1995 ) for Menzel's case B recombination theory. 



3.2. Theory of the CO bands 

The energy of a diatomic molecule in rotational level J and 
vibrational level v can be expanded in the following way (Dun- 
ham p932al Il932b|) 



E{v,J) = he 



k,l 



k.l 



(6) 
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The parameters Yk.i are for the CO molecule taken from Far- 
renqet al. (1991). The first overtone bands result from coupled 



vib-rot transitions in the ground electronic state and obey the 
selection rules Av = 2 and A J = ±1. 

In the following, we assume the CO gas to be in LTE. Then, 
the levels are populated following a Boltzmann distribution. 



N,,j = -{2J+l)e^^ 



(7) 



where N and T are the total number density and temperature of 
CO molecules and Z is the total partition function, given as the 
product of the vibrational and the rotational partition function 



Z = Zy ■ Zj = J2 



(8) 



To account for optical depth effects, we calculate the line 
intensities from the transfer equation 

4 =B,(T)(l-e---). (9) 

The optical depth is = J n^, ds with the absorption coeffi- 
cient per cm 

f2j + l iV„ 



2j' + 1 N,j 



- 1 ^coM (10) 



where $co('^) is the profile function of the CO gas. The Ein- 
stein coefficients, Ayj-y ji, are from Chandra, Maheshwari & 
Sharma ( |1996| ). 




2.2925 



2.293 



2.2935 
A [/im] 



2.294 



Fig. 5. Individual vib-rot lines that form the 2-^0 band head. 
The solid curve is the (2, 51) (0, 50) transition. The dotted 
curves are the neighbouring transitions with rotational quantum 
numbers decreasing in steps of one; for the dashed curves, the 
quantum numbers increase in steps of one. 

Let Aiiiin denote the onset of the (2 0) bandhead. We de- 
fine it to be the wavelength where the intensity of the (2,51) 



(0, 50) line (see Fig.g) has dropped by a factor a:; ~ 5 from its 
maximum at Ao; the exact number is not important. We deter- 
mine Amin from Figure |to (2.29285 ± 0.00005) ^m. 
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Fig. 6. The onset of the spectrum above noise is choosen as the 
minimum wavelength Amin- 

The value of Amin depends of course on the various mech- 
anisms of line broadening. Most processes lead to a gaussian 
or nearly gaussian profile, like thermal (turbulent) motion of 
mean velocity vth (wturb), detector resolution Urcs, spherical 
wind Wwind (see e.g. Hamann & Simon 1986). Their super- 
position results again in a gaussian line of width Wgauss = 



"wind + "res- 



4 + ^^t'urb ■ 

Disk rotation, on the other hand, produces a double-peaked 
profile. In a gas ring rotating at velocity Wjot a line with lab- 
oratory frequency I'o is shifted up to a maximum frequency 
i^o • (1 + Wrot/c). In a mixture of disk rotation and gaussian 
broadening one finds the following expression for the minimum 
wavelength. 



1 



An 



1 



1 



(11) 



where visr = 8 km/s is the velocity shift of MWC 349 with re- 
spect to the local standard of rest. This relation only holds if 
the emission in the rise of the spectrum is optically thin. The 
thermal velocity vth of CO molecules at 3500 K is less than 
~ 1 km/s and negligible in comparison with the spectral reso- 
lution tires = 10 . . . 15 km/s. The same can be said about the 
turbulent velocity which is not likely to exceed the speed of 
sound. 

Eq. (11) can now be used to relate Uwind to the rotational 
velocity Wjot- The result is shown in Figure ^ The dashed lines 
above and below the solid line are for the lower and upper limit 
of Amin- If the CO gas performs no rotation, the gas must em- 
anate from the disk in a wind with about ~50 km/s, the same 
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A^i, = (2.29285 ± 0.00005) /im 
■ 1 5 km/ s 
Visr = 8 km/s 




20 40 
Vrot [km/s] 

Fig. 7. The dependence of Uwind on v^ot after Eq. ( plj ) is shown, 
a; = 10, and the remaining parameters are chosen as indicated 
in the Figure. The solid Hne is for Amin = 2.29285 /xm, the 
dashed line above (below) is for the lower (upper) limit of Amin- 



velocity as for the Pfund lines. Without a wind, one derives 
a rotational velocity of ~ 58 ± 7 km/s. This value is larger 
than other rotational velocities around MWC 349 observed up 
to now, the highest being 42.5 km/s for He I (Hamann & Simon 



1986) 



4. Modeling tlie CO band emission 

4.1. A gaussian line profile for the CO gas 

First, we neglect rotation and assume a pure gaussian line pro- 
file. The velocity of the CO gas is taken from Figure^ as 
I'gauss « 50 km/s. 

For the calculation of the optical depth we assume the 
absorption coefficient to be constant along the line of sight. 
Then is simply given as the product of the absorption co- 
efficient per CO molecule, n^, times the CO column density, 

T^ = K^- Nco- 

In Fig.|^ we plotted for temperatures between 3000 K 
and 5000 K, and a gaussian velocity of the gas of 50 km/s. In 
the wavelength range around the 3 — > 1 band head (2.322- 
2.325 pm) the absorption coefficient is higher than in the region 
around the 2^0 band head (2.293-2.295 fj,m). This effect be- 
comes stronger with increasing temperature. Thus, the turnover 
from optically thin to optically thick emission depends not only 
on the column density but also on the temperature: the 3^1 
band head becomes optically thick at lower temperature and 
lower column density than the 2^0 band head. 

Next, we vary the CO column density and fix Tco 31 
4000 K. The synthetic spectra in Figure ^ are calculated for the 
case of extreme optical depth (upper line) and for CO column 
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Fig. 8. Absorption coefficient per CO molecule, k^, for tem- 
peratures of 3000 K (dotted hne), 4000 K (dashed line), and 
5000 K (solid line). The gaussian velocity of the CO gas is 
50 km/s. 
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Fig. 9. Calculated CO band emission spectra for a temperature 
of the CO gas of 4000 K, a gaussian velocity of ^ 50 km/s, and 
a CO column density increasing from bottom to top: 10^", 5 • 
10^", 10^^, and 2 • 10^^ cm^^. The upper line represents the 
case for r oo. 



densities (from top to bottom) of 2 • 10^^ 10^\ 5 • 10^", and 
10^" cm^^. With increasing column density the optical depth 
enlarges and the intensity tends towards its limiting blackbody 
value. In addition, the spectra show two characteristic features: 
a varying intensity ratio of the two band heads and a broaden- 
ing of the band head structure. 
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Fig. 10. The flux density ratio of the 3^1 and 2^0 band 
heads as a function of temperature and column density. The 
ratio is highest in the opticaUy thin case (A'^co < lO'^" cm^^) 
and approaches ~ 0.96 for r ^ oo. 



1 - 



E 
a. 



0.5 




2.293 



2.294 



A [/in 



Fig. 11. Theoretical (2 0) CO band head emission for a 
gaussian line profile with Ugauss ~ 50km/s. The CO column 
density is 2 • 10^^ cm~^ (dashed line), 5 • 10^° cm^^ (solid Une) 
and 10^" cm~2 (dotted line). 



The first property can better be seen in Figure 10 where 
we plotted the flux density ratio of the 3 — > 1 and 2 — > band 
head with temperature for different column densities. This ratio 
decreases with increasing column density, i.e. increasing r, and 
approaches 0.96 for r — > cx). 



0.5 




2.295 



2.305 



Fig. 12. Comparison between theoretical (2 0) band head 
spectra of optically thin (upper panel) and marginally optically 
thick emission (lower panel). In both cases, the temperature is 
3500 K. 



The determination of the column density with only the help 
of this ratio is, however, not definite because the flux density 
ratio also depends on temperature. From our observations we 
get no information at all on the column density in this way be- 
cause the observed strength of the 3 — > 1 band head is subject 
to some uncertainty as discussed in Section ^. 

The second characteristic feature is the broadening of the 
band head structure with increasing column density, especially 
at the long wavelength side. In Figure [IT] the observed 2 
band head region (histogram) is overlaid on the theoretical 
spectra with column densities of 2 • 10^^ (dashed line), 5 • 10^" 
2.295 (solid line), and 10^" cm^^ (optically thin case, dotted Une). 

At first glance, it seems that the observations might also be 
fitted by optically thin emission, but the upper panel of Fig- 
ure |l2| which includes a somewhat larger wavelength interval 
is clearly inconsistent with the optically thin case whereas the 
observations can be reconciled with a synthetic spectrum of 
column density A^co — 5 • 10^° cm^^ (lower panel of Figure 
|l2| ). An additional disagreement is that the flux density ratio of 
the bandheads is highest for the optical thin case which can not 
be confirmed by the observations. 

So we take A^co — 5 • 10^" cm^^ as the best fit CO column 
density. The synthetic spectra of the 2-^0 band head look very 
similar for temperatures in the range of 3000 — 5000 K. For a 
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Fig. 13. Modelled CO band emission (black line) for a temper- 
ature of Tco — 3500 K and a gaussian velocity of ^ 50km/s 
overlaid on the observations (grey histogram). The model pa- 
rameters are given in Table |^. The insert shows a blow-up of 
the 2^0 band head. Lower panel: residue of observations 
minus model. Some yet unidentified emission lines seem to be 
present. 



2.314/im, whereas the features at 2.3253 and 2.3084 /im seem 
to be yet unidentified emission lines. 

For a distance to MWC 349 of 1.2kpc we find a CO emis- 
sion area projected on the sky of Aqo — 1-36 • 10^® cm^ and a 
mass of Mqo — 3.2 • 10^'' g. Our Aqo is about 5 times larger 
than the value of 2.5 • 10^^ cm^ found by Geballe & Persson 



(1987) who could only derive a lower limit. 



4.2. CO bands from a Keplerian rotating disk 

We investigate whether the CO bands in MWC 349 can alter- 
natively be explained by emission from the circumstellar disk. 
For demonstration, we first calculate the optically thin emis- 
sion from an infinitesimally narrow ring. For Kepler rotation, 
the orbital velocity is 



«™t = yjGAh/r. (12) 
and the line of sight velocity, vy^s, of a gas element in the ring 



Wrot • COS a ■ COS I 



(13) 



Here G denotes the gravitational constant, r the ring radius, 
A/* the stellar mass, i the incUnation angle of the ring, and 9 
the azimuthal angle. The resulting profile of the CO band head 
smoothed to the observational resolution of 15km/s shows a 
characteristic shoulder at short and a maximum at long wave- 
lengths relative to an unbroadened band head (see Fig. The 
distance between shoulder and maximum decreases with de- 
creasing orbital velocity. The CO profile of a rotating ring is 
evidently very different from what is observed in MWC 349. 



more sensitive limitation of the temperature range we must take 
a look at the 3-^1 band head. From there we conclude that 
Tco = 3500^000 K. Unfortunately, an exact determination of 
the CO temperature fails because the higher bandheads are not 
available within our high resolution spectrum. 

Table 3. Parameters for CO used for the calculations shown 
in Fig. |3) A}^^ denotes the foreground extinction, Aco and 
Mco are the emitting CO area projected to the sky and the 
mass of CO derived from the fit, respectively. 



A^co 


Tco I'gauss 




Aco 


Mco 


[cm-2] 


[K] [km/s] 


[mag] 


[cm^] 


[g] 


5 • 10^" 


3500 ~ 50 


10 


1.36 ■ 10^'' 


3.2 ■ 10^* 



With the derived values for column density and temperature 
we model the synthetic CO band spectrum taking into account 
the foreground extinction, Aj,^'^, and fitted it to the observa- 
tions (Figure [T^; Table p| summarizes the parameters used). 
The lower panel of Fig. O shows the residue of observed mi- 
nus calculated spectrum. Some remaining features may result 
from observational uncertainties (indicated as dashed parts), 
for example the dip around 2.3225 /im and the region around 



'in 





1 1 1 1 


1 j 1 1 1 1 j 1 1 1 1 1 1 1 

A Tco = 3500 K ^ 
1 v,„, : km/s - ■- ' 

I \ 20 km/s 

1 40 km/s ■ 

1 / \\ 60 km/s 






1/ )K\ "J-intr = 1 5 km/s 

■i / -^a' 




^ / 

/ / 






/ / ;■ 






/ / /■ 




/ / 






/ / ■■' 











2.293 2.2935 2.294 2.2945 

A [/im] 



Fig. 14. Optically thin CO band emission from rings of con- 
stant rotational velocity (Vrot — 60, 40, 20 and km/s) and for 
an inclination angle of i = 0°. The spectra are smoothed to the 
observational resolution of 15 km/s. 
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We next calculate the emission from a whole disk viewed 
nearly edge-on, using the radiative transfer of Eq. (^. We as- 
sume power laws for the temperature and surface density, 



T{r)=T,-[- 



and 



= Tn • - 



S(r) = S 



• 



= Sn- 



ap 



2q 



(14) 



(15) 



r.i is the inner radius of the Kepler disk where the rotational 
velocity has its maximum value Vi. We fix the inner radius 
by putting Vi — 60km/s (see Figure^. Further, we assume 
To = 5000 K, which is the dissociation limit of CO molecules, 
and an inchnation angle of 10°. The radius 7-out, out to which 
the CO (2^1) band is thermally excited, is a free parameter; 
the velocity there equals Vout- For these disk calculations we 
surmize that the CO bands at 2 micron are for radii greater 
than Tout only subthermally excited due to lower densities. 
Therefore their contribution to the overall spectrum is proba- 
bly small. The CO disk is geometrically flat so that along each 
line of sight through the inclined disk the density and tempera- 
ture are constant. The total flux from the disk is obtained from 
a straight forward numerical integration. 



Table 4. Parameters for the fit of the CO disk calculations 
shown in Fig.|l5[ Fixed parameters are the mass of 26A/q, 
To = 5000 K, p = 0.5, 9 = 1.5, and Vi = 60km/s which leads 
tor,- ~ 6.4 AU. 



"out 


Tout 


T(rout) 


iVo 


iV(rout) 


[km/s] 


[AU] 


[K] 


[cm-2] 


[cm-2] 


40 


14.4 


3333 


8.5 ■ 10'" 


2.5 ■ 10'" 


30 


25.6 


2500 


5.8- 10'" 


7.3 ■ 10'^ 



By a systematic variation of the free parameters 
(So, Wout,P, ?), we find that the shape of the spectrum is 
strongly influenced by the column density at the inner edge and 
by the size of the disk, whereas the exponents p and q play a 



minor role; we therefore fix them following Hayashi (1981) to 
p = 0.5 and q — 1.5. S imilar numbers are derived for the dust 
disk (Ki-aus et al. |2000|). 



The disk models then contain only two free parameters: the 
outer radius, rout, and the column density at the inner edge, 
A'o. To explain the observed band head profile with a rotat- 
ing disk, one needs a large spread in the rotational velocity be- 
cause otherwise one cannot surpress the shoulders seen in Fig. 
p^ . The large velocity spread implies a large emitting area, and 
therefore a low column density, so that the emission becomes 
very optically thin (r ~ 10""'), even at the inner edge of the 
CO disk. The outer velocity Wout must lie in the range between 
30 and 40 km/s. If Wout were higher, the shoulder of the band 
head would become too broad, if it were smaller, the profile 
would get too narrow. Our best models of a rotating disk are 



E 
4 



1 .0 



0.5 




Vj = 60 km/s 
V;, = 40 km/s 
= 30 km/s 



2.293 



2.294 
A [/im] 



2.295 



Fig. 15. Best fit to the observations (histogram) by Kepler disk. 
Model parameters from Table ^ 

presented in Fig. |l5[ Considering the quality of the data, they 
are marginally acceptable. 

We add that if one assumes dust and gas to be thermally 
decoupled, it is in these models nevertheless possible that the 
strength of the CO band head amounts to / ~ 5 % of the 
continuum, as observed, although the CO emission is very 
optically thin. As the dust emits like a blackbody, the equa- 
tion f B^{Td) = Tyi?iy(Tco) is approximately fulfilled for 
Td ~ 1000 K, Tco - 4000 K and r - 10"^. 



5. Discussion 

Mathematically, there is also the possibility to combine the 
wind model of Fig. 1 1 (solid line), which is based on a gaus- 
sian profile of ~ 50 km/s half width, with a Kepler rotation 
as long as v^ot < 40 km/s. This is shown in Figure |6[ where 
a rotational component with Vrot = 40 km/s (black curve) is 
added to the gaussian velocity profile. Otherwise the basic pa- 
rameters are the same as for the wind (LTE; CO column density 
- 5 • 102"cm-2;Tco = 3500. ..4000 K). 

Although the fit is again very good, this scenario meets the 
following difficulties: If the rotational velocity v^ot < 40 km/s, 
the stellar mass of 26 M© implies a distance greater than 
14 AU. If the CO gas is coupled to the dust, it should have the 
same temperature of ~ 950 K as the dust grains at that dis- 



tance (Kraus et al. 2000), which is much too low. If CO is lo- 
cated above the disk and thermally decoupled from the dust, 
it might achieve the required temperature of ~ 3500 K by ra- 
diative heating. But the CO column density of 5 • 10^° cm~^ 
then implies a gas surface density which is at least an order 
of magnitude greater than what follows from the dust disk 
(adopting standard conversion factors Sdust : Sh — 10~^, 
Nco ■■ ~ 10-4). 
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rotation 
axis 




observable 
ring sector 



Fig. 17. Sketch of the disk around MWC 349 seen edge on. Because of the small inclination angle, only emission from sector 1 
where the radial velocity is small will reach the observer Light from other places, for example, from the points marked 2, 3, and 
4 is absorbed by the bulge. 



1.0 



0.5 




2.293 



2.2935 2.294 
A [fim] 



2.2945 



Fig. 16. Best fits to the observed 2^0 band head. The grey 
line represents the pure gaussian fit of Figure 13 and the black 
line shows the combination of a rotational component and a 
gaussian velocity component whose velocities are determined 
according to Eq.[rT|. 



The scenario of a pure wind is very vague in its details, but 
gives a remarkably good fit. We might think of dense molecular 
blobs emanating from the disk. The total emission area of the 
blobs is ~ 0.6AU2. 

Table 5. Angles 9 and velocity ranges Urot.ios for the calculated 
sectors shown in Figure |l8|. The inclination angle was taken as 
10° and the rotational velocity is 98 km/s. 





en 


't'rotjlos 


[km/s] 


dashed line 


180 


-96.5 


..96.5 


solid line 


75 


-59.0 


..59.0 


dotted line 


33 


-27.6 


. . 27.6 



But there is still another alternative with A^co ~ 5 • 
10^" cm^^, which we discuss now. If the CO band emission 
arises in the disk, it must come from inside the evaporation 
radius of the dust, rovap, because of its high temperature of 
~ 3500 K. According to the Kraus et al. model of the dust 
disk, the distance is then smaller than 2.88 AU according to a 
rotational velocity greater than 90 km/s. Let us assume ?;rot — 
100 km/s. To avoid the problems with the line profile in case 
of such a high v^ot, as discussed in subsection 4.2 and demon- 
strated in Figure |T^ we propose a configuration as depicted in 
Figure |l7| The disk has a bulge and because of its small incli- 
nation angle i ^ 10°, this bulge blocks the light from the inner 
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2.293 2.294 2.295 

A [fJ,m] 



Fig. 18. Fit to the observed 2^0 band head for ring segments 
with a Keplerian velocity of v^ot — lOOkm/s. The geometry is 
depicted in Fig. 17. The angles and the corresponding velocity 
ranges are given in Table ^ The best fit is shown by the solid 
line. 



(see Fig. 13). The band head has in this case an optical depth of 
order unity (iVco — 5 • 10^° cm"^). 

Alternatively, the CO might form in a thin layer above the 
dust disk at radii between 6.5 and 15 AU. In this scenario, the 
CO emission is optically very thin (A^co — 6 • lO^^cm"^; 
T ~ 10^"*), and the gas is thermally decoupled from the dust. 
As we could only produce fits of mediocre quality (see Fig. 15), 
due to the characteristic line profiles under Keplerian rotation 
(see Figure this is not our favorite configuration. 

Finally, we suggest that the CO bands come from the inner 
edge of the circumstellar disk. Because of the required high 
temperature and column density, the CO gas must be located 
inside the evaporation radius of the dust, i.e. at distances less 
than ^ 3 AU from the star Although this would imply a high 
rotational velocity component which would usually lead to a 
characteristic shoulder in the band profile, which is not ob- 
served, the disk may well have a bulge, as depicted in Figure 
|l7[ This bulge absorbs the CO emission from the near part of 
the rotating inner disk edge, so that the observer sees only a 
sector on the far side where the radial velocities are all smaller 
than ^ 60km/s. With such a geometrical configuration, a sat- 
isfactory fit is also possible (see Figure |l8|). 

Acknowledgements. We would like to thank Frank Shu, the referee, 
for his helpful comments 



edge of the disk on the near side to the observer. We therefore 
receive emission only from sector 1 in Figure |7| where the ra- 
dial velocity is low. In Figure [l^ we show computations of the 
flux received from sector 1 for different opening angles 9. The 
ranges of radial velocities in the sectors are summarized in Ta- 
ble | 

6. Conclusions 

We present low and high resolution spectra in the near infrared 
of the B[e]-star MWC 349. The wavelength interval 2.285- 
2.342 /j,m contains mainly emission of the first overtone bands 
of the CO molecule as well as the Pfund series of atomic hydro- 
gen. From modeling the Pfund lines under the assumption that 
they are optically thin, we find that they come from the inner 
part of the H II region around MWC 349 and are in LTE. The 
hydrogen line profiles are approximately gaussian and thek 
width indicates a wind velocity of ^ 50km/s. 

The CO molecules are at a temperature of 3500^000 K, 
the population of the vibrational levels is close to LTE. The 
width of the 2^0 band head indicates a velocity broaden- 
ing of the order of 50-60 km/s, depending on the broadening 
mechanism, wind or rotation. The location of the hot CO gas 
producing the band emission could not be clarified, but we were 
able to propose several scenarios. 

Because the CO as well as the Pfund lines can be fitted with 
a gaussian velocity component of 50 km/s, the first possibil- 
ity is that the CO bands arise in the windy transition layer be- 
tween the disk and the Hll region, presumably in dense clumps 
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